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Abstract 

The idea of the functional renormahzation group and one-loop improved 
renormalization group flows are reviewed. The associated flow equations and 
nonperturbative approximations schemes for its solutions are discussed. These 
techniques are then applied to the strong interaction in the framework of an 
effective quark meson model which is introduced in great detail. The renormal- 
ization group analysis of the two flavor quark meson model is extended to finite 
temperature and quark chemical potential which allows for an analysis of the 
chiral phase diagram beyond the mean field approximation. 



1 Introduction 

Quantum Chromodynamics (QCD) is the quantum field theory of strongly in- 
teracting matter. Many of its vacuum features have been tested experimentally 
over a wide range of momentum scales. At high momenta or small distances 
the asymptotic freedom of non-Abelian gauge theories can be used to apply 
perturbative methods for the computation of physical observables. Due to the 
running of the strong gauge coupling the situation becomes significantly more 
complicated at smaller momentum scales or larger distances. Here perturbation 
theory breaks down and nonperturbative methods are called for. There are no 
analytical methods starting from first principles which allow to treat QCD at 
larger distances, where the strong gauge coupling becomes large and pertur- 
bation theory fails. The main reason for this difficulty is that QCD describes 
qualitatively different physics at different length or energy scales. 

The situation is further complicated at finite temperature and/or baryon 
density. For instance, at very high temperatures perturbative calculations are 
plagued by serious infrared divergences. Furthermore, it is generally expected 
that at high enough temperature and densities hadronic matter attains a state 
in which chiral symmetry is almost restored and its fundamental degrees of 
freedom, the quarks and gluons, are no longer confined. The system undergoes 
a phase transition from the ordinary hadronic phase to a chirally restored and 
deconfined quark giuon plasma (QGP). Furthermore, recent theoretical studies 
reveal an increasing richness in the structure of the phase diagram of strongly 
interacting matter. 

Therefore nonperturbative methods are indispensable to obtain quantita- 
tively reliable results because one has to deal with large couplings and a pertur- 
bative expansion of physical observables fails. One such nonperturbative method 
is given by the renormalization group (RG) . The RG method represents an effi- 
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cient way to describe critical phenomena and phase transitions. It can be used 
to characterize universal and non-universal aspects of second-order as well as 
first-order phase transitions and is very well adapted to reveal the full phase 
diagram of strongly interacting matter. In the context of the phase diagram 
the RG-method has been applied in [H [21 [3l [4] using a two-flavor quark-meson 
model, which captures essential chiral aspects of QCD. 

In the following we will introduce the concept of the effective average action 
and its associated renormalization group equation. Since it is clearly impossible 
to review here the many facets of the renormalization group developments we 
will concentrate on the so called proper-time renormalization group approach 

mm- 

2 Renormalization group methods 

The renormalization group deals with the effect of a scale change in a theory. 
The central issue is the understanding of the macroscopic physics at large dis- 
tances or at low momenta in terms of the underlying fundamental microscopic 
interaction. In order to understand the evolution from the microscopic to the 
macroscopic scales one has to consider the quantum or statistical fluctuations 
on all scales in between. The general RG idea is to treat the fluctuations not 
all at once but successively from scale to scale [7 . 

This RG idea combined with functional methods yields the so called func- 
tional RG. By means of functional methods the computation of generating func- 
tional of correlation functions becomes feasible. All important physical infor- 
mation is contained in the correlation functions once the fluctuations have been 
integrated out. Instead of evaluating correlation functions by averaging over all 
fluctuations at once, only the change of the correlation functions induced by 
an infinitesimal momentum shell of fluctuations is considered. This goes along 
Wilson's philosophy of integrating out modes momentum shell by momentum 
shell. 

From a technical point of view this means one has to work with functional 
differential equations - the so called RG or flow equations - instead of func- 
tional path integrals which is usually the case in standard quantum field theory. 
The differential structure of the RG equations has a larger versatility and of- 
fers several advantages compared to an integral formulation. It is analytically 
and numerically better accessible and more stable. This is of great interest in 
non-Abelian gauge theories such as QCD. During the evolution from the mi- 
croscopic to the macroscopic scales these theories turn from a weak to a strong 
coupling becoming thus nonperturbative at macroscopic scales. In this sense RG 
methods provide an analytical and powerful tool to investigate nonperturbative 
phenomena of quantum field theories and statistical physics. 

There are various RG methods known in the literature [H El [TOl [Til lEl [Ml 
[14] . One particular formulation of RG flows is based on the concept of the 
effective average action Ffe, which is a simple generalization of the standard 
effective action F, the generating functional of the one-particle irreducible (IPI) 
Green functions 15J. The generalization is achieved by implementing an infrared 
cutoff scale k in the functional integral that defines the effective action. F^ is 
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then obtained by integrating over all modes of the quantum fields with Euclidean 
momenta larger than the infrared cutoff scale, i.e. > k^. In the limit A: — > 0, 
the infrared cutoff is removed and the effective average action becomes the full 
quantum effective action F containing all fluctuations. For any finite infrared 
cutoff k the integration of quantum fiuctuations is only partially done. The 
influence of modes with momenta < k^ is not included. 

On the other hand, in the limit k ^ oo (or to some finite ultraviolet cutoff A) 
the effective average action matches the bare or classical action, which does not 
contain any fluctuations. Hence, the knowledge of the fc-dependent F^. allows 
to interpolate from the bare action in the ultraviolet to the full effective action 
in the infrared. As the scale k is lowered more and more quantum fluctuations 
are taken into account. As a consequence F^ can be viewed as a microscope 
with varying resolution length scale ^ 1/fc. It averages the pertinent fields over 
a d-dimensional volume with size l/k"^ and permits to explore the system on 
larger and larger length scales. This is similar to a block-spin transformation 
on the sites of a coarse lattice where more and more spin-blocks are averaged 
over. 

The dependence of the effective average action F^ on the scale k is governed 
by the flow equation 



Here, F^ ' [0] stands for the second functional derivative with respect to some 
given field 4>. The trace involves an c?-dimensional integration over momenta 
(or coordinates) as well as a summation over internal indices (e.g. fiavor, color 
and/or Dirac indices). It is convenient to introduce the logarithmic variable 
t = ln(A:/A) which induces the equivalence of the partial derivatives dt = kd/dk. 

The infrared cutoff function which regularizes the generating functional is 
denoted by Rk and depends on the IR scale k. This regulator term, defined in 
momentum space, can be viewed as a momentum-dependent mass term. The 
regulator function is not completely arbitrary but has to satisfy certain condi- 
tions in order to provide an IR and UV regularization for the flow equation. 
Since it screens the IR modes in a mass-like fashion it should always be a posi- 
tive function regarding its arguments. Furthermore, it has to vanish for A: — > 
removing the regulator from the flow. This condition ensures that the full effec- 
tive action is recovered in this limit. Finally, it has to diverge for k ^ oo (or for 
A: — > A) . This limit ensures that the correct initial condition lim Tk = S\ in the 



ultraviolet is reached where denotes some bare action. This last condition 
acts as a delta functional in the defining integral representation of F^ which is 
dominated by the stationary point of the action. This limit justifies the use of 
a saddle-point approximation which filters out the classical field configurations 
and the bare action. 

The flow equation ([T]) is an exact equation in the sense that it can be derived 
from first principles and is often labeled as Exact RG (ERG). It has a simple 
one-loop structure but is not of one-loop order because it depends on the full 
field-dependent inverse average propagator F^ [(pj+Rk. This one-loop structure 




(1) 



k^A 
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has some advantages because technical comphcations due to overlapping loop 
integrations do not arise. 

The solution of the flow equation corresponds to a trajectory in the theory 
space, the space of all possible action functional which are spanned by all pos- 
sible invariant field operators. The starting point of the trajectory in this space 
is the bare action and the end point the full effective action F. As already 
mentioned, the explicit form of the regulator function Rk is not fixed uniquely. 
Different choices for Rk will lead to different trajectories thus reflecting the RG 
scheme dependence of the flow. Of course, the end point of all trajectories in the 
theory space is independent of the choice of the explicit form of the regulator 

Rk- 

General methods for the solution of the functional flow equation are rare. 
Nevertheless, one can turn the flow equation into a coupled system of non-linear 
ordinary partial differential equations for infinitely many couplings by expanding 
the effective average action in terms of invariants. In order to reduce the infinite 
system to a numerically manageable size one needs to truncate the most general 
form of Ffe. This can be accomplished in various systematic approximation 
schemes. One such approximation scheme is the 'operator expansion' which 
constructs the effective action from operators with increasing mass dimension. 

An example for this type of truncation is the gradient expansion where 
is expanded in numbers of derivatives. For a scalar field theory with one real 
field one obtains 

FfcM = j d'x + \Zkmd^^f + . (2) 

The first (constant) term Vk corresponds to the scalar effective potential and the 
first correction includes the field-dependent wave function renormalization Z^. 
In general, it is not guaranteed that the chosen truncation scheme converges. 
One has to control the truncation error separately as for any expansion. 

The flow of the effective average action, Eq. ([T]), is not the only exact one. 
It is closely related to other well-known exact flows like e.g. Wegner-Houghton 
flows [Hj or Callan-Symanzik flows [T71 HH] . 

Another class of RG flow equations can be obtained by approximations to 
Exact RG flows. An example of such an approximation is the so-called Proper- 
time RG (PTRG) flow which we want to investigate in the following. Compared 
to the ERG flow the proper-time flow has a numerically simpler and physically 
more intuitive representation and yields in the lowest-order gradient expansion 
smooth and analytical threshold functions. 

2.1 Proper-time RG flow 

Proper-time flows in their original formulation are one-loop improved RG flows 
[3 HI E] • They can be obtained by an RG improvement in an one- loop flow 
equation [HI 1101 HH [12] . Since the PTRG flows do not depend hnearly on the 
full field-dependent propagator they cannot be exact flows [23] . Furthermore 
one can show that the (non-exact) PTRG flow starts to deviate from standard 
perturbation theory already at two-loop order |24j . 
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Nevertheless, within a background field formalism a generalized PTRG flow 
can be formulated which is exact. This is achieved by adding further terms to the 
standard PTRG flow J5j. In the following an ad hoc derivation of the standard 
one-loop improved PTRG flow is presented. The discussion is restricted to an 
one-component scalar theory with a general interaction. The generalization to 
other theories is straightforward. 

2.1.1 Derivation of the standard PTRG &ow 

The starting point for the derivation of an one-loop improved RG flow is the 
one-loop effective action 

r]^'°°^m - s,,,.M + IttIuS^J^Im , (3) 

where the trace denotes a sum over all momenta. The classical action is la- 
beled by S'ciass[</'] and its second fleld-derivative by •S'^fassI'^] ~ '^^5'ciass[0]/<50<5(?!). 
The ill-defined logarithm is regularized by means of a Schwinger proper-time 
representation 

oo 

ri-'°°P[<A] = 5e,a..M -I J ^/fc(r)Trexp(-r41[^]) , (4) 



where the regulator function /a;(t) provides an IR cutoff k. The lower (upper) 
limit of the proper-time integral has inverse mass dimension two and regularizes 
the UV (IR) region of the theory. Since the momenta and the proper-time 
variable r are coupled, high momenta correspond to small proper-time values 
and vice versa. Taking the fc-derivative of Eq. (jlj the standard PTRG flow 
equation is obtained as 

oo 

dtTki^] ^at/fe(r)Trexp(-Tri^)[0]) , (5) 



where the classical action S^^!^^^ on the right hand side has been replaced by the 

scale-dependent effective action F^^-* . This represents the RG improvement [H] . 

Similar to the ERG flow the regulator function fk has to satisfy some con- 
ditions: Since the flow should start from the initial action Fa in the UV it is 
required that lim /fe(r) = 0. The second condition lim fk=io{T) — regularizes 

fc— »A T^OO 

the IR since the upper proper-time limit coincides with lower momenta. Finally, 
at the end of the evolution, the condition lim /fe(r) = 1 reduces the proper-time 

k — *0 

regularization to the usual Schwinger proper-time regularization. 

One class of regulator functions fl^\T) with ti > 0, which fulfill all these 
required conditions can be expanded in terms of incomplete F-functions and 
have the structure 
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The fc^-term in the exponential acts as a mass term and controls the IR behavior. 
Inserting this basis set of regulators in Eq. ^ and performing the proper-time 
integration yields the flow 

C\ n+l 
^5r^- — I ■ (^) 

For n > the PTRG flow cannot be an exact flow since it does not depend 
linearly on the full propagator. For n = the flow depends linearly on the full 
propagator and one could conjecture that it is an exact flow. In fact, it has the 
form of a Callan-Symanzik flow, but it is not exact. Furthermore, the flow looks 
like an ERG flow with a mass insertion Rk — fc^ but is not precisely an ERG 
flow. In contrast to the ERG flow, the momentum integration is not regularized 
in the UV because the condition lim Rk ^ oo is not fulfilled in this case. All 

momenta contribute to the flow and an additional UV rcnormalization would be 
required in contrast to the ERG flow. In addition, for this flow the Wilsonian 
single momentum-shell interpretation becomes questionable. 

Nevertheless, the standard PTRG flow is a well-defined approximation to a 
first-principles flow like the ERG flow It is an approximation to a back- 

(2) 

ground field ERG fiow where additional terms proportional to dtT^, are ne- 
glected and the difference between fluctuation flelds and background fields has 
been omitted. With these additional terms the generalized PTRG flow becomes 
exact. So far, applications of the PTRG flow are typically based on flows of 
the form of Eq. ([7]) and further approximations thereof. The standard PTRG 
flow with the lowest-order gradient expansion truncation has a simple polyno- 
mial momentum dependence which will lead to simple flow equations. However, 

(2) 

it is still an open issue how these additional flow terms proportional to SfF). 
which make the flow exact affect the results obtained with the standard PTRG 
flow. First results in this respect for an 0{N) symmetric scalar theory indicate 
a minor influence of these addition flow terms at least at criticality in three 
dimensions [211 HZ] . 

In the following the standard PTRG flow will be used to investigate strongly 
interacting matter which is described by QCD. 

2.2 Applications 

In QCD the quarks and gluons represent the microscopic degrees of freedom, 
whereas the macroscopic degrees of freedom are the observed color neutral parti- 
cles like the mesons, baryons and/or glueballs. Hence there must be a transition 
from the microscopic to the macroscopic degrees of freedom. Therefore, in QCD 
the relevant degrees of freedom change with the scale k which is suitable for an 
RG treatment. In order to apply RG techniques to QCD an initial starting 
effective action has to be formulated. 

When constructing effective models for these macroscopic degrees of freedom 
one usually relies on the guiding symmetries of QCD because a first-principle 
derivation from QCD is still missing. One important symmetry of QCD is the 
local SU (Nc) color invariance which is related to confinement. This symmetry 
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cannot be used here since the observed hadronic spectrum consists of color bUnd 
states. We will concentrate on the chiral dynamics of QCD and consider QCD 
with only two light quark flavors. To a good approximation the masses of these 
two flavors are small compared to the other quark flavors. For vanishing current 
quark masses the classical QCD Lagrangian does not couple left- and right- 
handed quarks. Omitting the axial anomaly and baryon number conservation 
it exhibits a global chiral invariance under the SU^iNj) x SUii{Nf ) symmetry 
group where Nf denotes the number of quark flavors. In the observed hadron 
spectrum only the vector- like subgroup SUviNf) is realized which implies a 
spontaneous symmetry breaking of the chiral SUL{Nf) x SUbX^/) symmetry 
down to the SUv{Nf) symmetry. This symmetry breaking predicts for Nf = 2 
the existence of three light parity-odd Goldstone bosons, the three pions tt. 
Their comparably small masses are a consequence of the explicit chiral symmetry 
breaking due to the finite current quark masses. 

At scales above ~ 2 GeV the dynamics of the relevant degrees of freedom 
appear to be well described by perturbative QCD. At somewhat lower scales 
quark condensates, bound states of quarks and gluons emerge and confinement 
sets in. 

To each such non-perturbative phenomenon one can associate an appropriate 
scale. Focusing on the physics of scalar and pseudoscalar mesons and assuming 
that all other bound states are integrated out it appears that all these scales 
are rather well separated from each other. The compositeness scale k^, where 
mesonic bound states are formed due to the increasing strength of the strong 
interaction is somewhere below 1 GeV. The chiral symmetry breaking scale 
at which the chiral quark condensate develops a non-vanishing value will be 
below the compositeness scale, typically around 500 MeV. The last scale where 
confinement sets in, is related to the Landau pole in the perturbative evolution 
of the strong coupling constant and is of the order of Aqcd ^ 200 MeV. 

For scales between the ranges k^ < k < k^ the most relevant degrees of 
freedom are quarks and mesons and their dynamics in this regime is dominated 
by the strong Yukawa coupling g between them. This picture legitimates the 
use of a quark meson model only if one assumes that the dominant QCD effects 
are included in the meson physics. Below the scale k^^ the strong coupling as 
increases further and quark degrees of freedom will confine. Getting closer to 
Aqcd it is not justified to neglect those QCD effects which certainly go beyond 
the meson dynamics. Of course, gluonic interactions are expected to be crucial 
for an understanding of the confinement phenomenon. But due to the increase 
of the constituent quark masses towards the IR the quarks decoupled from the 
further evolution of the mesonic degrees of freedom. As long as one is only 
interested in the dynamics of the mesons one expects that the confinement 
on the mesonic evolution has only little influence even for scales below Aqcd- 
Hence there are good prospects that the meson physics can be described by an 
effective quark meson model [28l EH |30] . 

The Ansatz for the effective action at the compositeness scale fc^ is given by 
the Nf = 2 quark meson model 

Ffe^fe, = / d^x g{a + i-f^rn)) q + ^{d^a)^ + {d^^f - f/(a,fr)| (8) 
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where the purely mesonic potential is defined as 

The quark fields are denoted by q and g is the Yukawa coupling which de- 
scribes the interaction strength of the isoscalar-scalar a and the three isovector- 
pseudoscalar pion fields tt with the quarks. The value for g is usually fixed to 
the value of the constituent quark masses Mq which are of the order of 300 
MeV in the vacuum. Here we have neglected the S'C/(2)-isospin violation and 
consider averaged quark masses. The explicit symmetry breaking parameter c 
in the mesonic potential controls the value of the pion mass. For c = the 
chiral limit is obtained and the pion mass vanish. In the chiral limit the action 
is invariant under global chiral SU{2)l x SU{2)fi symmetry transformations. 
The last two parameters of the model, A and , are fixed to the mass of the 
sigma field and the pion decay constant 

This effective action at the compositeness scale emerges from short dis- 
tance QCD in basically two steps. Starting from QCD in the UV one first 
computes an effective action which only involves quarks. This corresponds to 
an integration over the gluonic degrees of freedom in a quenched approximation. 
This will generate many effective nonlocal four and higher quark vertices and a 
nontrivial momentum dependence of the quark propagator. In the second step 
decreasing the scale further these four and higher quark interactions will cause 
the formation of mesonic bound states. Thus, at the compositeness scale not 
only quarks but also composite fields, the mesons, are present and interact with 
each other. The four quark interactions have been replaced by mesonic fields. It 
is obvious that for scales below the compositeness scale a description of strongly 
interacting matter in terms of quark degrees of freedom alone would be rather 
inefficient. 

The Ansatz for the effective action ([8]) incorporates a truncation to four 
quark interactions. Higher quark interactions have been suppressed at the com- 
positeness scale. This leads to a purely quadratic mesonic potential at with 
a positive mass term. This implies that the chiral symmetry is restored at the 
compositeness scale and the scalar expectation value {a) will vanish 

The effective action also corresponds to a gradient expansion to lowest order 
(cf. Eq. ^ where wave-function renormalizations are neglected. It serves as 
the required initial action for the RG evolution. The evolution is started at the 
compositeness scale where the mesons are assumed to be formed due to strong 
gluonic interactions. During the evolution towards the IR this effective action 
for the quark meson model is used as a truncation for an approximate solution of 
the RG flow. Instead of expanding the mesonic potential in a polynomial form 
for the fields we will be more general here and allow for arbitrary higher mesonic 
O (4) -symmetric self-interactions in the potential. We will solve numerically the 
flow equations for the full effective potential. But the running of the Yukawa 
coupling will be neglected. 

The dynamics in the beginning of the scale evolution just below the compos- 
iteness scale is almost entirely driven by quark fluctuations. These fluctuations 
rapidly drive the squared scalar mass term in the action to negative values. 
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This then immediately leads to a potential minimum away from the origin such 
that the vacuum expectation value (cr) becomes finite. This happens at the 
chiral symmetry breaking scale k-^ < k^, not far below k^. The reason for this 
behavior lies in the suppression of the meson contributions. All meson masses 
are much larger than the constituent quark masses around these scales and are 
therefore further suppressed during the evolution. Furthermore, the quarks are 
strongly coupled to the mesons since the Yukawa coupling is relatively large. 
Below fc^ the systems stays in the regime with spontaneous chiral symmetry 
breaking. Around scales of the order of the pion mass the evolution of the po- 
tential minimum stops. The reason for this stability of the vacuum expectation 
value is that the quarks acquire a relatively large constituent mass Mq. These 
heavy modes will decouple from the further evolution once the scale drops below 
Mq. The evolution is then essentially driven by the massless Goldstone bosons 
in the chiral limit. Of course, for non- vanishing pion masses the evolution of the 
model is effectively stopped around scales k ~ ttItt. Quarks below such scales 
appear to be no longer too important for the further evolution of the mesonic 
system. Because of confinement effects quarks should anyhow no longer be in- 
cluded for scales below Aqcd- The final goal of such an evolution is to extract 
phenomenological quantities like meson masses, decay constants etc. in the IR. 
This can be achieved in a straightforward computation from the effective action 
r = lim Ffc. 

fc-»0 

So far we have considered the vacuum fluctuations that contribute to the flow 
of the effective average action F^. Using certain phenomenological quantities in 
the IR as input we have flxed all model parameters at the compositeness scale. 
The extension of the presented RG formalism to a system in thermal equilibrium 
system with finite net-baryon or net-quark number density is straightforward. 
In such systems the effective average action plays the role of the grand canonical 
potential Q and depends also on the temperature T and on one averaged quark 
chemical potential fiq if again isospin symmetry is assumed. 

In the Imaginary-time or Matsubara formalism a finite temperature T re- 
sults in (anti-)periodic boundary conditions for (fermionic) bosonic fields in the 
compact Euclidean time direction with radius 1/T. This leads to a replacement 
of the zeroth component of the momentum integration with discrete (even) odd 
Matsubara frequencies for (fermions) bosons. Thus, for finite temperature a 
four-dimensional quantum field theory can be interpreted as a three-dimensional 
system plus an infinite tower of Matsubara modes for each degree of freedom. 
Since all Matsubara frequencies are proportional to the temperature all massive 
Matsubara modes will decoupled from the dynamics of the system at high tem- 
perature. One therefore expects an effective three-dimensional theory with the 
bosonic zero modes as the only relevant degrees of freedom. This leads to the 
phenomenon of dimensional reduction. 

The implementation of a finite quark chemical potential in the quark me- 
son model is also straightforward. One has to add a term proportional to 
i^q J d'^xqjoq to the Ansatz ([HJ for the effective average action. This Ansatz 
together with the choice of a half- integer regulator function /^,"\ see Eq. 
where the variable n is replaced by 3/2 is used as the input to derive the proper- 
time fiow equations. The choice of a half-integer power in the regulator function 
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corresponds to a three-dimensional momentum cutoff and allows for an analytic 
evaluation of the Matsubara sums. This is one advantage of using the proper- 
time flow with this regulator compared to other RG flows. Finally, one obtains 
for the scale-dependent grand canonical potential the proper-time flow equation 



with the pion- E.„ = .JW+2VtF^, the cr-meson = y/W^^2^V^7+H^I, 
and quark energies Eq = ^JW^+g^^ . The primed potential denotes the (fy^- 
derivative of the potential, i.e., ilj, :— dilk/dcj)^ and correspondingly the higher 
derivatives. The potential for this model depends on the expectation value 
of the square of the chiral 4-component field 0^ which coincides with (cr)^ 
since (tt ) ~ 0. The scale-dependent effective meson masses are defined as 

J, = 2il'^. + 40^f2J! and rn^ f. = 20!^. where the potential has to evaluated 
at the global scale-dependent minimum (jy^ = 0q. The dynamically generated 
constituent quark mass Mq = g(j)Q is proportional to the minimum since the 
running of the Yukawa coupling is not considered here. 

The summations of the fermionic (bosonic) Matsubara sums yield analytical 
functions which represent the corresponding mass threshold functions for finite 
temperature and quark chemical potential. The threshold functions are an 
important nonperturbative ingredient generally appearing in flow equations. 
They control the smooth decoupling of massive modes from the evolution once 
the IR cutoff scale k drops below the corresponding mass. The different degrees 
of freedom contribute in an additive way to the flow. One recognizes the three 
degenerate pion, the sigma and the quark/ antiquark threshold functions in the 
square brackets. The flow equation has an overall scale factor fc^ which reflects 
the correct dimension of the effective potential in d = 4 dimensions. This can 
be seen explicitly by rewriting all mass threshold functions in a dimensionless 
form. The fermionic contributions enter with a negative sign due to the fermion 
loop and have a degeneracy factor of (2s + l)NcNf with s = 1/2. The quark 
chemical potential enters only in the quark/antiquark threshold functions with 
the appropriate sign as it should be. It influences the bosonic part of the flow 
equation only implicitly through the meson masses. The flow for the vacuum can 
be deduced immediately of the full flow © by examining the limits T = 0, = 
of the finite temperature and quark chemical potential threshold functions. In 
this limit all hyperbolic functions tend to one and only the vacuum threshold 
functions remain. 

The flow equation ([9]) constitutes a coupled, highly non-linear, partial dif- 
ferential equation which can be integrated numerically in principle in two ways: 
either one discretizes the unknown potential f2fc on a (/)^-grid or expands the po- 
tential in powers of 0^ around its minimum 0q . The advantage of the potential 
expansion is that only a finite set of coupled flow equations has to be solved, 
depending on the chosen expansion order. In particular, this yields the beta 
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functions for the couplings of the potential. For each higher order of the poten- 
tial expansion, however, a new coupled beta function is introduced increasing 
the numerical effort drastically. A further drawback is that the potential is only 
known around the minimum 0q once the system has been solved [22 [H] . This 
is different for the grid solution: Here, the potential is not only known around 
the minimum but also for arbitrary c/)^. This is of importance, for example, 
in a first-order phase transition where two degenerate minima of the potential 
emerge. In this case the knowledge of all local minima is required to describe the 
phase transition correctly. This is cumbersome in a potential expansion, except 
for some simple potentials. With an additional explicit symmetry breaking term 
in the potential every potential minimum has always a finite value because the 
symmetry is never restored exactly. A precise determination of the critical tem- 
perature of a first-order transition is very difficult within an expansion scheme 
around only one potential minimum. 

Since we have opted for the grid solution, the field cfP' has been discretized 
for a general potential term on a regular grid. Thus, for each grid point a 
flow equation is obtained which leads finally to a coupled closed system. As 
initial condition at the chiral symmetry breaking scale we use a symmetric 
potential Ansatz. All model parameters are fixed in such a way that they match 
the physical IR vacuum quantities accordingly. 

For zero chemical potential and for two massless quark fiavors a second-order 
phase transition is found in which the spontaneously broken chiral symmetry is 
restored. The phase transition and fixed point structure belongs to the 0(4)- 
universality class with the corresponding critical exponents. The value for the 
critical temperature is a non-universal quantity and depends on the model pa- 
rameters. For finite pion masses, i.e. with an explicit symmetry breaking term 
in the potential, the second-order phase transition is washed out and turns into 
a smooth crossover. This results in a shift of the "pseudo-critical" temperature 
defined as the inflection point of the order parameter towards larger values. 

For flnite chemical potential, the 0(4)-universal second-order phase transi- 
tion persists up to a tricritical point which is a critical point where three phases 
coexist. The curvature of the second-order transition line Tc(/i) has a negative 
slope. Thus, the location of the tricritical point must be below the value of the 
critical temperature at zero chemical potential. The tricritical point belongs to 
a trivial Gaussian flxed point with mean-fleld critical exponents. The precise 
location of this tricritical point is in general not known and again depends on 
the model parameters [31]. Thus, the existence of this point, the shape of the 
transition line and its universality class are predictions within the underlying 
quark meson model. At higher chemical potential and smaller temperatures 
the phase transition changes initially to a single flrst-order phase transition. 
For smaller temperatures we observe a splitting of the transition line and two 
phase transitions emerge. The left transition line represents a first-order tran- 
sition down to the T — Q axis. At this transition the order parameter jumps 
not to zero but to a finite value. The chiral symmetry remains spontaneously 
broken and is only restored for higher chemical potentials which then produce 
the second (right) transition line. At this right transition line wc initially find 
a second first-order transition where the order parameter jumps to zero and 
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chiral symmetry is restored. But for smaller temperatures close the /z-axis the 
order parameter tends smoothly to zero. It seems that this transition is again 
of second-order. If this is true we infer that there must be a second tricritical 
point in the phase diagram in the chiral limit. 

As the pions become massive the tricritical point turns into a critical end 
point (CEP). For temperatures below the CEP a first-order curved transition 
line is found which persists down to the /x-axis. For each value of the explicit 
symmetry breaking parameter (finite pion mass) , there is a corresponding CEP. 
In an extended, three-dimensional (T,/x,to^) phase diagram these points arrange 
into a critical line. The static critical behavior of this line falls into the uni- 
versality class of the Ising model in three dimensions corresponding to the one- 
component scalar 0'*-theory in three dimensions. As the pions become massive 
no symmetry remains which would require the order parameter to have more 
than one component. On the other hand, in the chiral limit we have a chiral 
SU{2)l X SU{2)r symmetry restoration which is isomorphic to the 0(4) sym- 
metry. In this case the order parameter should be made up of four components, 
one sigma and three pions. 

In the chiral limit below the splitting point, the right second-order transition 
turns into a crossover for finite quark masses. Analogously, the second tricritical 
point, if it exists, should turn into a critical point. Some remnants of this critical 
point can indeed be seen in the vacuum expectation value and in the behavior 
of the meson masses. But a detailed analysis of this point is postponed to a 
future work. 

3 Summary 

In this lecture note we have presented an brief introduction to the functional 
renormalization group methods with a focus on the flow equation for the effective 
average action. We specialized and discussed in great detail the standard proper- 
time flow and its relation to the effective average action. In the second part we 
used this proper-time RG to explore the QCD phase diagram within an linear 
quark meson model for two quark flavors. This model captures essential features 
of QCD such as the spontaneous breaking of chiral symmetry in the vacuum 
and can therefore yield valuable insight into the critical behavior associated with 
chiral symmetry. 
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